Large next-to-leading order QCD corrections to pentaquark sum rules 
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Using configuration space techniques, we calculate next-to-leading order (NLO) five-loop QCD cor- 
rections to the correlators of interpolating pentaquark currents in the limit of massless quarks. We 
obtain very large NLO corrections to the spectral density which makes a standard sum rule analysis 
problematic. However, the NLO corrections to the correlator in configuration space are reasonable. 
We discuss the implications of our results for the phenomenological sum rule analysis of pentaquark 
states. 
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INTRODUCTION 

There is a continuing interest in exotic states of strong 
interactions that differ from the standard nonexotic 
mesons and baryons (see e.g. 0, 0)- Nonexotic three- 
quark baryons have been intensively studied in |^ 
Finite mass effects were studied in 5, 6, 7] where the 
next-to-leading order (NLO) QCD corrections to the cor- 
relators of finite mass baryonic currents were determined. 
This led to an improved precision of the sum rule analy- 
sis. Different aspects of the physics of multiquark states 
in QCD have been discussed long ago It is notewor- 
thy that the analysis of multiquark exotics can provide 
valuable information on the details of quark-gluon inter- 
action relevant to nuclear physics 0, flOl . 

Candidates to be described in QCD as multiquark 
states are the deuteron 0, ^| or Jaffe's dihyperon 
H 0, Since the dihyperon state can only decay 

weakly, it would be expected to be quite narrow. To the 
best of our knowledge the dihyperon state is the first state 
to attract attention in the modern context of QCD. The 
discovery of gluon bound states would be a triumphant 
confirmation of QCD and would allow for a quantitative 
check of QCD in the completely new sector of the glue- 
ball states Another class of multiquark states, 
the pentaquark states, have recently become the focus of 
intense theoretical and experimental studies. 

The properties of all of these multiquark (or 
multigluon) states can be studied in a model independent 
way through the method of the QCD sum rule analysis. 
In the QCD sum rule analysis one analyzes the operator 
product expansion of current-current correlators of inter- 
polating local fields which have the quantum numbers of 
the multiquark states under study. It is well known that 
radiative corrections have a strong impact on the results 
of the sum rule analysis. In this paper we derive the nec- 
essary tools that allow one to compute the as radiative 
corrections to multiquark sum rules in the limit when 
the quarks (or antiquarks) are massless. As a specific 
example we apply our method to pentaquark correlators 
and calculate the NLO radiative corrections to the pen- 
taquark current correlator and the spectral density of a 



specific pentaquark current. 



CALCULATION 

An important ingredient in the formulation of the op- 
erator product expansion (OPE) analysis for the pen- 
taquark states is the choice of the interpolating current. 
The result depends strongly on the choice of the interpo- 
lating current as has been pointed out before in the sum 
rule analysis of the dibaryon 8] . The same holds true for 
the sum rule analysis of pentaquark states 
detailed analysis of the dependence on the interpolating 
current is out of the scope of a short letter, it will be pub- 
lished elsewhere '2^ . In the following we shall present the 
tools needed to calculate the NLO QCD corrections for 
interpolating quark currents of any composition. 



Generalities 

In the QCD sum rule analysis the prime object of study 
is the correlation function 



(1) 



where the interpolating current j{x) is a local operator 
with the quantum numbers of the pentaquark baryon 
state O. It has to be constructed from four quark and 
one antiquark fields such that its projection onto the pen- 
taquark state |0(p)} is nonzero: 



(OU(O)ie(p)) = Ae, P^^ml 



(2) 



Since the interpolating current j{x) is not unique, one 
immediately faces the question of the optimal choice for 
the interpolating currents. Recall that the problem of 
choosing the optimal current already arose in the case of 
baryons 0, 01 where the currents are constructed from 
three quark fields. In the pentaquark case there are 
five quark (antiquark) fields to build the interpolating 
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currents and correspondingly the number of indepen- 
dent currents with the correct quantum numbers is much 
larger (see 0| and references therein). 

The treatment of the interpolating current j{x) in its 
most general form, i.e. in the form of linear combinations 
of all possible independent local operators with the cor- 
rect pentaquark quantum numbers, is a very unwieldy 
problem. In this pilot study of the NLO corrections to 
interpolating pentaquark currents we limit ourselves to 
the study of only one simple interpolating currents with 
the required pentaquark quantum numbers. As in the 
case of mesons and baryons, we construct our interpolat- 
ing current j(x) from quark fields without derivatives. 

Tools for the NLO calculation 



The LO diagram for the pentaquark correlator is given 
by the four- loop diagram FIG.^a). There are two types 




FIG. 1: LO contribution (a) and examples for the NLO prop- 
agator (b) and dipropagator corrections (c) 



This choice of renormalization scale avoids the appear- 
ance of ln(47r) and je factors in configuration space cal- 
culations. The relation of fix and the usual renormaliza- 
tion scale /i of the MS-scheme is given by 

The dipropagator two-loop amplitude for a pair of quarks 
with open Dirac indices leads to an integral that is a bit 
more difficult to calculate. The result for the dipropaga- 
tor correction including the LO term reads 

5'2(a;)|NLO = Sa{x^)^{Y'x^,(g>Yxu 
a 1 
47r s 

+ a3rf^«)r3 a/3%a:.)} (5) 
where the coefficients oi, bi and 03 are given by 

11 , 1 11 

ai = -l-—e, bi = -l--e, ag = -- - -e, 

and where 

^liri'^i" -YTin- (6) 

Eqs. (O and jSJ allow one to calculate the NLO cor- 
rections to n-quark(antiquark) current correlators of any 
composition using purely algebraically algorithms with- 
out having to compute any integrals. 

Renormalization of the interpolating current 



of NLO five- loop corrections. First there are the one- 
loop corrections to single quark propagators S{x), one 
of which is shown in FIG. H^b). We shall refer to these 
corrections as the propagator corrections. Then there 
are the dipropagator corrections connecting two different 
quark propagators, one of which is shown in FIC^c). 

It turns out that it is very convenient to calculate these 
corrections in configuration space, in particular if the 
quarks are treated as massless i21j. For the propagator 
correction Si{x) we obtain 



Siix) 



NLO 



S,{x)\^o[l~CF^\{f?x^'Y 
So{x')rx,\^l-CE^\{fiWyY^ 



where in the Euclidean domain one has 

-iT{2 - e) 



27r2-^(a;2)2-e- 



(4) 



fix is a renormalization scale in dimensional regulariza- 
tion appropriate for calculations in configuration space. 
The space-time dimension is parametrized by I? = A — 2e. 



In a NLO calculation one has to account for mixing 
effects between operators when going through the renor- 
malization program. Mixing can occur when gluons are 
exchanged between the lines in the pentaquark correla- 
tor (dipropagator corrections). In order to keep track of 
fiavours we first treat the case of an interpolating cur- 
rent composed of five massless quark(antiquark) fields 
with different flavours. 
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where C is the charge conjugation matrix. The in- 
terpolating current consists of a baryonic part B = 
^ijkf^qiT (j^r^qk ^ mesonic part M = {(f^qi,)- If the 
gluon is exchanged within the mesonic part, the renor- 
malization factor is the usual one for the mesonic opera- 
tor (see e.g. [i^V 



Zm = 1 - 



7re 



(8) 



If the gluon is exchanged within the baryonic part, the 
renormalization factor is given by the known renormal- 
ization factor for the baryonic operator [j^ . 
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as 



(9) 
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These two contributions will be referred to as factoriz- 
ing contributions. When the gluon is exchanged between 
the mesonic and baryonic part one obtains new operators 
which do not have the structure of the original current. 
There are two contributions of this kind which we refer 
to as mixed contributions. If the gluon is exchanged be- 
tween the mesonic part and the [qf^Cq!^) piece of the 
baryonic part, we obtain a contribution of the operator 



Oi = e^^^[{q'^Ca^Pqi) + {q'^ Ca'^^ q\)]qt{q{a^P qi) 



where <7°'^ 



(10) 



|[7",7^]. If the gluon is exchanged be- 
tween the mesonic part and the g| piece of the baryonic 
part one obtains a contribution of the operator 

O2 = e^'\q^;'Cqi)a^P[q[{q{a^pqt) - \qWi^o.pql)]- 

(11) 

The renormalized current reads 

= ZMZsf'^'' - + ■ (12) 



NLO RESULTS 

In order to obtain NLO results, let us consider the 
interpolating current 



Adding Eqs. H15() and (|16|l one obtains the renormalized 
correlator 



TlRix^) = n(a;^) + An = 360 1 + — (3 + SL^^) 
-6 (1 + ^ (7 + UL,)) - 4 • 6 (1 + ^ (6 + 9i, 



(17) 



with = ln(a;2^3f ) = \n{x'^e^^'^ /A). 

The spectral density is obtained by calculating the dis- 
continuity of 



j = e«''=(u'^^Cd^)d"(s'=u'=) 



(13) 



io Jl-s(px)(x2)--^3-2s^^ (18) 

for the cases a — 7 and a — 7— e, where J\{z) is the Bessel 
function of the first kind lU]. One obtains (s — p^) 



. ^ , , ^ , oit f 617 

p{s)lpn{s) = 360( 1 + ^ + 3L 



with L = In(^Vs) and p^{s) = sV378000(47r)^ 



with the quantum numbers of the pentaquark baryon O. 
The current is constructed in such a way that it cannot 
directly dissociate into a neutron and a kaon since the 
respective color singlet parts have the wrong parity. 

The result of our NLO calculation for the bare corre- 
lator reads (j = e'''"'{d''C~'^v!''^)d''{u'' s^)) 



(0|rj(x)j(0)|0) = S^{x^f{x^f^^^x^IV{x^) (14) 



where 



n(x2) = 360 1 + —{ii\x 



-4.6|l + ^i(p5,.,T2)'(?-l| (15) 



The first line contains the factorizing NLO contibutions 
while the second and third lines contain the mixing con- 
tributions. The counter term contains contributions from 
the rcnormalization factors Zm and Zb as well as from 
the operators 0\ and O^- It reads 



cts \ 3 /II 28\ /9 
An = ^ <^ 360- - 6 ^|-4-6(--7 

TT £ \ e 3 



(16) 



DISCUSSION 

The first NLO contribution in both the correlator and 
the spectral density comes from the factorizing part of 
the diagrams, i.e. the case where baryon and meson cor- 
relators are multiplied in configuration space and do not 
mix by gluon exchange. This contribution dominates the 
final result but is probably irrelevant for physics as it 
does not contain a pentaquark bound state. The physi- 
cally relevant mixing contributions are smaller since they 
are suppressed both in the number of colours and a 
factor 4 coming from the evaluation of Dirac traces. 

The result for the spectral density ((T^ shows that the 
NLO corrections to the spectral density are large, where 
the NLO corrections to the factorizing parts are some- 
what smaller (note, however, that the relevant scale is 
not s but rather s/5 because of the number of lines in 
the correlator). This spoils the conventional momentum 
space QCD sum rule analysis. The NLO corrections to 
the correlator function (|17() in configuration space are 
more reasonable. This suggests a sum rule analysis in 
configuration space, based on the correlator l)17|l . How- 
ever, it is not clear whether the accuracy of such a sum 
rule analysis in configuration space will be sufficient for 
physical applications. 
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CONCLUSION 

We have calculated NLO perturbative QCD correc- 
tions for pentaquark current correlators which turn out 
to be very large. One has to conclude that the NLO QCD 
sum rule analysis of pentaquark states is fraught with dif- 
ficulties. This complicates the mass determination of the 
pentaquark states using a QCD sum rule analysis. One 
may have to take recourse to model calculations to de- 
termine the properties of pentaquark states such as the 
one based on chiral solitons [ijj . 
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